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Abstract 

We report density-matrix renormalization group calculations of spin gaps in the quarter-filled 
correlated two-leg rectangular ladder with bond-dimerization along the legs of the ladder. In the 
small rung-coupling region, dimerization along the leg bonds can lead to large enhancement of 
the spin gap. Electron-electron interactions further enhance the spin gap, which is nonzero for all 
values of the rung electron hopping and for arbitrarily small bond-dimerization. Very large spin 
gaps, as are found experimentally in quarter-filled band organic charge-transfer solids with coupled 
pairs of quasi-one-dimensional stacks, however, occur within the model only for large dimerization 
and rung electron hopping that are nearly equal to the hopping along the legs. Coexistence of 
charge order and spin gap is also possible within the model for not too large intersite Coulomb 
interaction. 

PACS numbers: 71.10.Fd, 71.30.+h, 71.27.+a, 75.10.Pq 
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I. INTRODUCTION 



Coupled pairs of one- dimensional (ID) chains, referred to as two-leg ladders, have been 
widely investigated within the spin 1/2 antiferromagnetic Heisenberg Hamiltonian because 
they exhibit behavior different from both the isolated ID chain as well as the two-dimensional 
(2D) lattice. Ignoring electron occupations, two different kinds of couplings between the ID 
chains, as shown in Figs. 1(a) and (b), have been of interest. We refer to the system of 
Fig. If rectangular ladder and that of Fig. 1(b) as a zigzag ladder. Within the spin 

1/2 Heisenberg Hamiltonian for the rectangular ladder, the energy gap between the lowest 
spin triplet (S=l) excitation and the spin singlet (S=0) ground state, hereafter the spin gap, 
is nonzero for arbitrarily small spin exchange along the rung bonds . *i 2 i 3 i 4 The zigzag ladder 
of Fig. 1(b) can also be considered as a ID chain with the interstack (intrastack) couplings 
corresponding to the nearest neighbor (next nearest neighbor) couplings in the ID chain. 
The spin gap in the antiferromagnetic zigzag ladder is nonzero only for the intrastack spin 
exchange above a critical valu o 5 ^ 1 18 . Nonzero spin gap is also found in the 1/2- filled band 
rectangular ladder within the Hubbard Hamiltonian with finite on-site electron-electron (e- 
e) interaction.^ Spin gaps and superconducting pair-pair correlation functions in weakly 
doped rectangular ladders have been investigated within Hubbard and t-J models^. 

The strong interest in the theory of undoped and weakly doped 1/2-filled band ladders 
owes its origin to possible connections with theories of high temperature superconductivity 
in the cuprate Whether or not significant spin gap occurs in ladder structures at other 
commensurate bandfillings is also an interesting question. The particular commensurate 
bandfilling that has been investigated so far is 1/4. Exact diagonalization studies of finite 
1 / 4-filled rectangular ladders with e-e and Holstein electron-phonon (e-p) interactions have 
demonstrated the coexistence of charge-order and spin-Peierls states^. The magnitude of the 
spin gap in the thermodynamic limit due to such a transition was not calculated. Density 
matrix renormalization group (DMRG) calculations for the rectangular ladder within an 
extended Hubbard Hamiltonian with repulsive on-site interaction U and nearest neighbor 
interaction V have shown that the checkerboard charge-ordering of Fig. 1(c) occurs only for 
V larger than a critical value V C (U)^. The calculated spin gaps in the thermodynamic limit, 
obtained from extrapolations, are nonzero for t± < t, where t± and t are the one-electron 
hopping parameters along the ladder rungs and legs, respectively, in agreement with weak 
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coupling renormalization group calculations^. The magnitudes of the extrapolated spin 
gaps are, however, tiny. Furthermore, the uncertainties in the extrapolated spin gaps are 
rather large in this parameter regime due to finite-size effects, and are comparable to the 
magnitudes of the spin gaps.^ii For example, for U — 8, V — and hopping parameter t± 
= 0.7 (all in units of t) the spin gap is ~ 0.03 ± 0.01, while for the same U and hopping 
parameter t± but V = 2, the spin gap is ~ 0.01 ± 0.005. The spin gap also goes to for 
hopping parameter t± > t, and consequently, significant spin gap is absent in the regular 
rectangular 1/4- filled band ladder for repulsive U > V . Significant spin gap does occur in the 
zigzag 1/4-filled band electron ladder with both e-e and e-p interactions. For t<i > 0.5858t s 
the ground state is a Bond-Charge-Density wave (BCDW) (see Fig. 1(d)) with coupled 
lattice distortion, charge order and very large spin gap.— 

Experimentally, large spin gaps have been observed at low temperatures in several re- 
cently discovered organic charge-transfer solids with strong pairwise couplings between quasi- 
1D stacks of organic molecules, and weak interpair coupling a 15 ! 16 ' 17 ! 18 . In all cases, there 
occurs first a metal-insulator transition at high temperature that is accompanied by leg 
bond-dimerization, following which there occurs an insulator-insulator transition which is 
accompanied by the opening of a spin gap. The temperatures T$g at which the spin gap 
transitions occur, ~ 70 K in the charge-transfer solid (DT-TTF)2[Au(mnt)2} 15: , and ~ 170 
K in (BDTFP)2PF6(PhCl)o.3^, are unusually high compared to the spin-Peierls transition 
temperatures Tsp ~ 15-20 K in the ID 1/4-filled band charge-transfer solid o 19 ^ 20 . 

The very large Tsg in the the coupled-stack systems suggest that the mechanism by 
which the gap opens is probably unrelated to the intrachain spin-Peierls transition, which 
in the 1/4-filled band is tetramerization, or dimerization of the dimerization.— Based on 
the experimentally observed high temperature dimerization, the ladder model of Fig. 1(d) 
has been suggested for these system o 15 i 16 i 17 i 18 . Within this picture, each dimer unit cell acts 
as a single site in an effective spin 1/2 rectangular ladder with significant spin gap (see 
below). The mapping from the dimerized 1/4-filled band electron model to the spin model, 
or the occurrence of spin gap in the original 1/4-filled band model have, however, not been 
demonstrated. 

Yet another 1/4-filled band ladder system in which relatively large spin gap (34 K) has 
been observed is a'-NaV^Os, which consists of V-0 ladders coupled through weak V-V 
bonds, with the V-0 layers separated by layers of Na + ions2i. Charge order involving the V 
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sites and spin gap occur at the same temperature in this system.^ Theoretical description 
based on the V-only 1/4- filled d-band extended Hubbard model predict the checkerboard 
charge-order pattern of Fig. 1(c). for the V-ladders . 22 ' 23 The calculated spin gap within the 
checkerboard charge-order phase for hopping parameter t± > t, as is appropriat e 22 ! 23 for 
a-NaX^Os, is however zero^i. The spin gap here most likely originates from the coupling 
between the 1/4- filled ladde ra 23 i 24 i 25 i 26 , which can lead to effective dimerization along the 
legs of the individual ladders with the checkerboard charge order (see below). 

In view of the interest in theoretical models for coupled 1/4-filled band two-chain sys- 
tems that give large spin gaps, we have performed DMRG calculations for the 1/4-filled 
band rectangular ladder with leg-bond dimerization within the extended Hubbard model. 
The purpose of our work is not necessarily to explain the observed experimental behavior 
of the organic charge-transfer solids or a'-NaX^Os, but to determine the plausibility of the 
application of the models to these systems. The minimum requirement of the applicability 
of the model is that significant spin gap occurs within it for realistic rung hopping parameter 
t±. We have considered two distinct parameter regimes, (i) V < Vc(U), corresponding to ho- 
mogenous charge distribution on the sites, and (ii) V > V C (U), which gives the checkerboard 
charge order in the infinite ladder.— In (i), we have considered mostly hopping parameter 
t± < t, as the results for larger t± are easily anticipated from these calculations. In (ii) 
we have considered hopping parameter t± > t because of the possible applicability of to 
a/-NaV205, in which the rung hopping is larger than the hopping along the leg bonds. In 
both cases we find spin gaps that are nearly an order of magnitude larger than those in the 
undimerized ladders. 

The paper is organized as follows. In section II we present the theoretical model and 
methodology. In section III we show the numerical results for the two cases, nearest neighbor 
electron-electron interaction V < Vc(U) and V > V C (U). Finally, in section IV we present 
our conclusions and discuss possible applications of the present work to the 1/4-filled band 
paired-stack charge-transfer solids and a'-NaX^Ov 

II. THEORETICAL MODEL AND METHODOLOGY 

We have performed DMRG calculations for the dimerized rectangular ladder of Fig. 1(d), 
within the Hamiltonian, 
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i i,A=-l,+l 

+V± ^ n^xUi -i 

i 

i,a,X=— 1,+1 

+t±X)(ct ljffCij _ v + /i.c.) (1) 

t,(T 

Here c- Acr creates an electron with spin a = |, j on site 2 in ladder leg A = ±1, 
n i,\,<r = c] A a Ci t \ tU , and 5 is the bond- alternation parameter along the legs. The parame- 
ters U, V, V±, t, t± have their usual meanings. In the following we present all quantities with 
dimensions of energy in units of hopping parameter t. We are interested in the parameter 
space U > V, V± > 0, and only in possible spin gap, as the conditions for charge gap and 
charge order have been discussed extensively in previous work 9,11 . For simplicity, we present 
results for V = V±_ and refer to both interactions as V. The spin gap A 5 is defined as usual as 
E(l) - E(0), where E(S) is the lowest energy in the state with total spin S. Note that for the 
1/4-filled band spin gap As = for hopping parameter t± = even with bond-alternation 
parameter 5^0. Although we are mostly interested in determining whether the spin gap is 
nonzero for arbitrarily small hopping parameter t±_ and bond-alternation parameter 5, our 
calculations span both t± < 1 and t±_ > 1. 

Before discussing our numerical results we present here brief discussions of the the physical 
mechanism behind the opening up of the spin gap in the dimerized ladder. Consider first 
the U, V,V± = band limit of Eq. d The electronic structure in this limit and for bond- 
alternation parameter 5 = is given by noninteracting bonding and antibonding bands 
split by 2tj_. For t±_ < 1 and bandfilling of 1/4, electrons occupy both bands and there 
are four Fermi points. For t± > 1, the lower band is 1/2-filled and the upper band empty. 
U > introduces a charge gap now, but the spin gap continues to be zero. Nonzero leg 
bond dimerization 5 ^ in the limit U, V, V± — opens a gap at the wavevector k = it /2a 
(where a = lattice spacing in the ladder leg direction). For small t± and bond-alternation 
parameter 5, both bands are again occupied. For £j_ ~ 1 now there appears a charge and 
spin gap 2(t± — 1 + 5) in the spectrum. For t± » 1, the energy gap at the Fermi surface 
in the bonding band is 45. These known results suggest that there exists a range of t± and 
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S where the spin gap can conceivably be significant for nonzero U, V, V±. 

The above physical picture does not explain the logic behind the proposed relation- 
ship between the dimerized 1/4-filled band model of Eq. (1) and the spin or 1/2-filled 
band ladders . 15 ! 16 ' 17 ' 18 This can be understood starting from a different limit, large bond- 
alternation parameter S ~ 1. In this limit each dimer unit forms a molecule with bonding 
and antibonding molecular orbitals (MOs) at energy ±2t. The single electron within each 
dimer molecule occupies the bonding MO, which forms an effective single site, as the unoc- 
cupied antibonding MO is too high in energy and does not play significant role in the low 
energy behavior of the system. The singly occupied bonding MOs, coupled through the weak 
intersite leg bonds and the rung bonds then constitute the effective sites of a 1/2-filled band 
ladder, for which we expect nonzero spin gap for arbitrarily small hopping parameter t±. 
Whether or not the spin gap is nonvanishing also for small 5 can be found only numerically. 
We will show that the spin gaps that we obtain numerically are considerably more enhanced 
compared to those obtained in the U = V = limit, viz., 2(t± — 1 + 5) for t± ~ 1 and 48 
for t± » 1. We will also show that the U-dependence of the spin gaps for the dimerized 
1/4-filled band ladder are very similar to those of the 1/2-filled band ladder. 

We calculated the spin gap for rectangular ladders up to 2 x 64 sites using the infinite- 
system DMRG algorithm with open boundary condition. Instead of a single site, a single 
rung is added to the building block each time. The number of states we keep is m = 300 for 
the S=0 state and 600 for the S=l state. We have confirmed that the accuracy of the singlet 
state is comparable to that of the triplet with these m. We checked some of our results 
against published results, 9 and found that the errors in our calculations are comparable, 
estimated to be less than a few percent. 

In principle, our calculations should be done with the finite-system DMRG algorithm, 
which is known to be more accurate than the infinite-system DMRG algorithm. The former 
procedure, however, requires considerably larger amount of time for each different ladder size 
and parameter set. The particular problem we consider requires that calculations are done 
with many different U, V, t± and bond-alternation parameter S (see below), with multiple 
system sizes L for effective L — > oo extapolations. The enormous amount of computational 
time that would be necessary for performing all the calculations that are reported below 
makes the finite-system DMRG algorithm impractical in the present case. We note, however, 
that in nearly all the cases we discuss below the calculated spin gaps are large and easy to 
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detect. We also note that the primary purpose of our work is to determine whether or not 
significant spin gaps are obtained within the dimerized ladder, and the precision required is 
not the same as within models where the gaps are tiny.— We have calculated a few of the 
spin gaps for several different system sizes, for representative parameters that cover both 
V < V C (U) and V > V C {U) within the finite-system DMRG algorithm using the same number 
of states as the infinite-system algorithm, and we compare these values with those obtained 
with the infinite-system algorithm in Table 1. As seen from the Table, the accuracies of the 
infinite-system algorithm are acceptable for our purpose, in both parameter regimes. 



III. NUMERICAL RESULTS 

A. Case 1. V < V C (U) 

In Fig. 2(a) we show our calculated spin gap As for 1/4-filled 2xL ladders, for the 
representative case of U = 8, V = 0, t± = 1 and < 5 < 0.05. For this moderately large 
the accuracy of the numerical results is very high and reliable results are obtained for 
the smallest bond-alternation parameter 5. The kink in the spin gap As plot for bond- 
alternation parameter 5 = 0.05 is real, - similar behavior has been seen previously for the 



undimerized 1/4-filled band ladder (see Fig. 3 in reference I 111). The dashed lines indicate 
the L — > oo extrapolations of the spin gaps, obtained by fitting the calculated spin gap 
As against a polynomial in 1/L up to 1/L 2 . For bond-alternation parameter 5 = 0.02 and 
0.05, we have also performed the extrapolations by retaining terms up to 1/L 3 . The two 
different extrapolations in each of these cases give the uncertainties in the spin gaps, which 
are 0.121±0.008 and 0.244±0.015 respectively. For all other cases fittings with 1/L alone 
gave extrapolations that are indistinguishable from the plots shown on the scale of the figure. 
The L —>■ oo extrapolated spin gap for bond- alternation parameter 5 = is zero in Fig. 2(a), 
in agreement with previous work.— For bond-alternation parameter S as small as 0.01, the 
extrapolated spin gap is nonzero for t± = 1. Furthermore, for this t± the spin gap at U 
= is 25. As seen in the Fig., for all bond-alternation parameter 5, the extrapolated spin 
gap As for U = 8 is larger. We therefore conclude that the spin gap is enhanced by the 
on-site e-e interaction, and is hence nonzero for the smallest bond-alternation parameter 5 
at t± = 1. The inset in Fig. 2(a) shows the weakly sublinear behavior of spin gap As against 
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bond-alternation parameter 5. 

Unfortunately, as the system size approaches the thermodynamic limit, the spin gap As 
decreases at smaller t± and the numerical accuracy of the DMRG results at large L are 
reduced. This has also been seen in previous calculations.^^ Even for the relatively simpler 
case of the 1/2-filled Hubbard ladder precise calculations of the spin gap are difficult for 
t± < 1.— Our calculations of the spin gap for the smaller t± = 0.8 in Fig. 2(b) are therefore 
for moderate bond-alternation parameter 5 > 0.05. Considering that at U = the spin 
gap As = 2(t± — 1 + 5) are zero for bond-alternation parameter 5 < 0.2, we see that the 
extrapolated spin gaps are once again strongly enhanced by U. Taken together with the 
weak-coupling result that the spin gap is nonzero (albeit small) for all t± < 1 at 5 = 0,— the 
results of Fig. 2(a) and (b) suggest that for nonzero 5 and U the spin gap is nonzero for all 
t±. The main difference between 5 = and 5 ^ is that the spin gap in the latter case is 
largest for t± > 1, exactly the region where the spin gap is zero for 5 = 0. The appearance 
of nonzero spin gap for the smallest bond-alternation parameter 5 indicates that there is 
indeed a qualitative similarity between the the 1/4-filled band dimerized rectangular ladder 
and an effective spin ladder. 

In Fig. 2(c) we have shown the effect of nonzero V on the spin gap for U = 8 and t± 
= 0.8 and bond-alternation 5 = 0.05. V enhances the spin gap very strongly, such that 
even for t± = 0.8 the spin gap is now quite large. The enhancement of the spin gap due 
to V can be understood physically, given the ^-dependence of spin gap Ag. Defining the 
bond-orders along the ladder legs in the usual manner, B it \ = (c-,a,o- c «+i,a,o- + h.c), we note 
that for nonzero bond-alternation 5 the absolute value of the difference between consecutive 
bond orders, AB = \B^x — -B,;+i,a| > 0. AB is then a wavef unction measure of leg bond- 
dimerization, 14 and for fixed bond- alternation 5 we expect the spin gap to increase with 
AB. We have calculated the exact AB(V) at fixed U = 8 for finite 1/4-filled periodic ID 
rings and ladders with t± = 0.0 and 0.8 with 12 and 16 sites. In Fig. 2(d) we have plotted 
AB{V) / AB(0) for U = 8 against V. As a consequence of the well known 4n vs. 4n+2- 
electron effect, we expect the curves for the infinite system to be bounded by the curves 
for the 12- and 16-site periodic systems for both tj_. The enhancement of AB(V) with V 
indicates enhanced effective leg bond-dimerization with V. Thus V increases the spin gap 
because it increases the dimerized nature of the ground state wavefunction. For U = 8, the 
checkerboard charge order of Fig. 1(c) is obtained^ for V > V C (U) ~ 2.6 at bond-alternation 



8 



parameter 5 = 0.— . Our calculations in Fig. 2(c) and (d) are for V < V C (U). We discuss 
V > V C (U) separately. 

In order to get more complete understanding of the behavior of the spin gap within 
Eq. Q w e have calculated the L — > oo extrapolated spin gaps for several sets of U, V, t± 
and 5. The parameter region of large t± > 1 and large 5 > 0.2 will clearly yield large 
spin gaps and is hence uninteresting. We are interested primarily in the t± < 1 region, 
thought to be applicable to the organic charge-transfer solids . 15 i 16 i 17 i 18 Also in the context 
of the charge-transfer solids, 5 < 0.1.— In Fig. 3(a) we have plotted the extrapolated spin 
gaps As versus t± for several different U, but V = 0, for bond- alternation parameter 5 = 
0.05. Because of the large uncertainties in the calculated As for small t± our results in 
Fig. 3(a) are limited to t± > 0.6. The results for the same U but larger bond- alternation 
5 = 0.1 are shown in Fig. 3(b). Fig. 3(c) shows the spin gap behavior for fixed U = 8, with 
three different V < V C (U). In all cases, and for all t± we find considerable enhancement of 
the spin gap by e-e interactions, as is summarized in Fig. 3(d). The behavior of the spin 
gap as a function of U in Fig. 3(d) is qualitatively similar to that of the uniform 1/2-filled 
band rectangular Hubbard ladder (see Fig. 4 in reference 0), suggesting again that the spin 
excitations of the dimerized 1/4-filled rectangular ladder can be obtained from those of the 
rectangular spin-ladder (note that increasing atomic U also increases the Coulomb repulsion 
between two electrons occupying the same bonding MO of the dimer unit^I. It has been 
claimed that such mapping can give rung spin exchange in the effective spin ladder larger 
than the nearest neighbor spin exchange along the legs, even for smaller hoppings along the 
rungs within the quarter-filled band electron ladder . 15 i 16 > 17 For the dimerized linear chain 
the transformation from the atomic basis to the dimer MO basis^ gives t/2 as the effective 
hopping between dimer MOs along the chain. This remains true here for the leg bonds. On 
the other hand, within the atomic basis of Eq. (1), there are two rung hoppings between the 
dimer units. The overall rung hopping between the dimer bonding MOs then continues to 
be t± after transformation to the dimer MO basis. Hence for all rung hopping t± > 0.5 in 
the 1 / 4-filled ladder the rung spin exchange in the effective spin ladder is larger than the leg 
spin exchange. We have confirmed this numerically: the L — > oo spin gap of the 1/4-filled 
band ladder at a fixed U with any specific t±/t is reproduced in the 1/2-filled band ladder 
with same U for much larger t±/t. 
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B. Case 2. V > V C (U) 



The calculations reported in the above are for V < V C (U) and t± < 1. For V > V C (U) 
and U,V » 1, the checkerboard charge order of Fig. 1(c) is obtained in the infinite ladder, 
while in finite systems there occur charge-charge correlations corresponding to this charge 
order.— The "occupied" sites in Fig. 1(c) behave as charge e particles with spin 1/2. At 
5 = 0, nonzero spin gap can occur only for small t±, where a spontaneous Majumdar-Ghosh 
dimerization spin gap 5 i 6 i 7 i 8 opens up.— There can, however, be no spin gap in this case for 
t_i_ > 1. The situation is different for nonzero 5, as indicated in our numerical results for 
large t± = 1 in Fig. 4 (a). V C (U) here for S = is 2.6^, and nonzero 5 is expected to slightly 
enhance V C (U). For t± = 1, we find then coexistence of the checkerboard charge order and 
spin gap in the region 3 < V < 4. With even larger t± = 2 the spin gap varies only weakly 
with V. The energy gap between the bonding and antibonding dimer MOs are now very 
large, and the system behaves here as a simple ID dimerized Heisenberg chain. 

The width of the region over which there occurs coexistence of charge order and spin gap 
clearly depends on the bond-alternation. Indeed, had we performed our calculations self- 
consistently with e-p couplings that modulated the hopping integrals along the leg bonds, 
as opposed to with rigid bond-alternation 5, the amplitude of the bond-dimerization would 
have decreased with V in the V > V C (U) region and the coexistence region would have 
been narrower. This is known from calculations within the ID 1/4- filled band e-p coupled 
extended Hubbard model where the bond dimerization decreases with V in the region V > 
V C (U)2&. The same physics can also be captured from calculations of the bond orders B(V). 
In Fig. 4(b) we show plots of AB(V)/AB(0) for the parameters of Fig. 4(a), obtained from 
exact diagonalizations of 12 and 16-site periodic ladders with 6 and 8 electrons, respectively. 
As in Fig. 2(d), the plots for the infinite system are expected to lie in between the curves for 
the 12 and 16-site ladders for each t±. Thus for the infinite ladder, the bond order difference 
is expected to first increase with V at t± = 1 up to V ~ 2, following which it is expected 
to decrease in the 2 < V < 4 region. Similarly, for t± = 2 the bond order difference is 
expected to be nearly flat in this region of V, with the decrease occurring at even larger 
V. In both cases, the qualitative behavior of the bond order differences and the spin gaps 
are the same, emphasizing the strong role of the leg bond dimerization. Similar relationship 
between the dimerization and the spin gap has also been observed in a recent calculation 
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on coupled ladders, where the effective dimerization is a consequence of the charge order on 
neighboring ladders. 26 

IV. DISCUSSIONS AND CONCLUSIONS 

In summary, the spin gap within the 1/4- filled band rectangular correlated-electron ladder 
is nonzero for the smallest bond-alternation, and increases monotonically with rung hopping 
t± in the region < t± < 1. The qualitative similarity between the dimerized 1/4- field band 
ladder and the spin ladder is easily understood in the large bond-alternation region, and 
persists at smaller bond-alternation, as is seen from the enhancement of the spin gap by 
dimerization, and the enhancement of the spin gap by both U and V for V < V C (U) at 
fixed dimerization. For V > V C (U), there can be coexistence of charge order and spin gap, 
especially for large t± > 1. 

In the context of the charge-transfer solids, the electronic parameters appropriate for 
quasi- ID cationic systems are known to be 2 ^ t ~ 0.1 eV, U/t ~ 6 — 8 and V ~ 2t < V C (U). 
Assuming similar intrastack parameters for the paired-stack charge-transfer solids, the very 
large Tsg i n these indicate spin gap As > 0.2. Assuming that the bond-dimerization 5 ~ 
0.1, as is true in most quasi-lD charge-transfer solids 2 ^, we see from Fig. 3 that such large 
spin gaps are obtained only for t±_ > 0.8, and certainly t± < 0.7 is impossible. Somewhat 
larger spin gap is obtained within the 1/4-filled band zigzag ladder for realistic e-e and e-p 
interaction parameters^, but the BCDW here appears only for td > 0.6. Thus both the 
electron ladder models would require substantial interchain hopping to give the observed 
large spin gaps. The difference between the two models is that while no charge order is 
expected within the rectangular ladder model for the charge-transfer solids, spin gap in 
the zigzag ladder is accompanied by strong charge order. Experiments that probe charge 
disproportionation are required to determine which of the two models apply to the real 
materials. 

As regards a / -NaV205, our numerical results of Fig. 4(a) are in agreement to the conclu- 



sions of references 



25 



and 



26l As shown in these papers, because of the coupling between 
the ladders in a'-NaV20s, the charge order can drive effective bond dimerization along the 
ladder legs that is manifested in the alternation of intraladder coupling parameters (there is 



a subtle difference between reference 



25l which models a'-NaV205 as a system with charge 
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order on only alternate ladders and reference l26l which assumes charge order in every ladder; 



this difference is not important in the present context.) The spin gaps in these references 
are explained within coupled dimerized spin chain models. DMRG calculations within the 
extended Hubbard model for two coupled ladders with charge order also find spin gap up 
on extrapolation^, but by necessity the maximum ladder length (20 rungs) is small here. 
Our calculations are thus complementary to previous work, showing that once it is assumed 
that interladder coupling gives rise to effective bond dimerization the essential physics is 
captured within the electronic Hamiltonian for the single ladder. 
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TABLE I: Spin gaps obtained by using the infinite-system and finite-system DMRG algorithms. 
In all cases U =8. 



Parameters 


L 


spin gap A s 






infinite 


finite 




40 


0.178337 


0.177095 


V=0 


48 


0.165737 


0.164359 


t ± =l, 5=0.02 


56 


0.157386 


0.155845 




64 


0.151652 


0.149933 




oo 


0.121±0.008 


0.111±0.006 




40 


0.235359 


0.234463 


V=2 


48 


0.234640 


0.233727 


t ± =0.8, (5=0.05 


56 


0.234241 


0.233316 




64 


0.234000 


0.233068 




oo 


0.233±0.001 


0.232±0.001 




40 


0.166847 


0.166470 


V=3 


48 


0.162380 


0.161969 


t ± =l, 5=0.05 


56 


0.159576 


0.159126 




64 


0.157700 


0.157209 




oo 


0.146±0.003 


0.145±0.003 
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Figure Captions 



Figure 1. Quarter- filled band two- leg electron ladders. Grey, black and white circles repre- 
sent mean electron occupations 0.5, > 0.5 and < 0.5, respectively, (a) Rectangular ladder 
with hopping integrals t and t± along the leg and the rung, respectively, (b) The BCDW 
broken symmetry in the 1/4-filled zigzag ladder, with average intrastack and interstack hop- 
ping integrals t s and td, respectively, (c) The checkerboard charge order in the 1/4-filled 
band rectangular ladder, (d) The bond-dimerized 1/4-filled band rectangular ladder. 

Figure 2 (a) and (b) Convergence behavior of spin gaps for the 1/4-filled band dimerized 
rectangular ladder for {7=8 and V=0. In (a) the rung hopping parameter t±=l, the maxi- 
mum number of rungs L max =80; O , A , □ and O correspond to bond- alternation parameter 
5=0, 0.01, 0.02 and 0.05, respectively. In the curves corresponding to 5=0.02 and 0.05, the 
upper L — ► oo extrapolation is using polynomial in ~ up to and the lower one up to j*. 

In (b), tj_=0.8; O , +and x correspond to (5=0.05, 0.1 and 0.2, respectively, (c) Spin gaps 
for {7=8, tj_=0.8, 5=0.05, with different V. O , A and ■ correspond to V=0, 1 and 2, respec- 
tively, (d) The normalized absolute difference in consecutive bond orders, AB(V)/AB(0), 
as a function of V for U = 8. Here O and □ correspond to periodic rings with 12 and 16 
sites (i.e., t±=0 in Eq. [T]and 5=0.05) , respectively; the • and ■ to ladders with 12 and 16 
sites (ij_=0.8, 5=0.1) respectively. In all cases the number of electrons is half the number of 
sites. 

Figure 3 (a) DMRG L — > oo spin gaps versus rung hopping t± for different U, for bond- 
alternation 5=0.05. (b) Same for 5=0.1. (c) Same for fixed {7=8 with different V. (d) 
Extrapolated spin gaps at 5=0.05 versus U, for three different t±. • and ■ correspond to 
V=l for rung hoppings ij_=0.8 and t±=l. 

Figure 4. (a) L — > oo spin gaps versus V for {7=8, 5=0.05 and two values of t±. For V > 3 
the spin gap coexists with the checkerboard charge order, (b) The normalized absolute 
difference in consecutive bond orders; The O and □ correspond to 12 and 16-site ladders, 
respectively, at t±=l. The • and ■ correspond to 12 and 16-site ladders, respectively, at 
t ± =2. 
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